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ABSTRACT 

Maps of the Galactic electron-positron annihilation radiation show evidence for 
three distinct and significant features: (1) a central bulge source, (2) emission in the 
Galactic plane, and (3) an enhancement of emission at positive latitudes above the 
Galactic Center. In this paper, we explore the possibility that Sgr A East, a very 
prominent radio structure surrounding the Galactic nucleus, may be a significant 
contributer to the central bulge feature. The motivation for doing so stems from a 
recently proposed link between this radio object and the EGRET 7-ray source 2EG 
J 1746-2852. If this association is correct, then Sgr A East is also expected to be a 
source of copious positron production. The results presented here show that indeed 
Sgr A East must have produced a numerically significant population of positrons, but 
also that most of them have not yet had sufficient time to thermalize and annihilate. 
As such, Sgr A East by itself does not appear to be the dominant current source 
of annihilation radiation, but it will be when the positrons have cooled sufficiently 
and they have become thermalized. This raises the interesting possibility that the 
bulge component may be due to the relics of earlier explosive events like the one that 
produced Sgr A East. 

Subject headings: acceleration of particles — cosmic rays — Galaxy: center — galaxies: 
nuclei — radiation mechanisms: nonthermal — supernova remnants 
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1. Introduction 

Electron-positron annihilation radiation from the Galactic Center region was first reported in 
the 1970s after the successful flight of a series of balloon-borne instruments (Johnson, Harnden 
& Haymes 1972; Johnson & Haymes 1973; Haymes et al. 1975). An unambiguous identification 
of the 511 keV line and the detection of the positronium continuum followed several years later 
(Leventhal, MacCallum & Stang 1978). Since that time, numerous follow-up experiments have 
provided additional monitoring of this source, though not always with a consistent measurement 
of the location or distribution of the emission, nor on the possible time variability of the 511 
keV line flux (see the review by Tueller 1993). However, the Oriented Scintillation Spectrometer 
Experiment (OSSE) on the Compton Gamma-Ray Observatory was recently able to carry out 
high-sensitivity observations of the Galactic plane and uniquely map the distribution of positron 
annihilation radiation and search for time variability of the emission (Purcell et al. 1997). 

The resulting maps show evidence for three distinct and significant features: (1) a central 
bulge source with a total flux (line plus continuum) of ~ 3.3 x 1CT 4 photons cm~ 2 s _1 and a 
Full Width at Half Maximum (FWHM) of ~ 4 degrees, (2) emission in the Galactic plane with a 
total flux of ~ 10~ 3 photons cm -2 s" 1 and a broad extension of more than 9 degrees FWHM in 
both latitude and longitude, and (3) an enhancement or extension of emission at positive latitudes 
above the Galactic Center with a total flux of ~ 9 x 10 -4 photons cm -2 s _1 and a FWHM of over 
~ 16 degrees. The average spectrum for the 1991 December viewing period of the Galactic Center 
(central bulge plus enhancement) exhibits a clear line feature at 509 ± 5 keV. The line width (2.5 
keV) is consistent with the instrumental resolution. The observed ratio of the 511 keV line flux to 
the total positronium flux corresponds to a positronium fraction of 0.98 ± 0.04. In addition, OSSE 
found no evidence for time variability. 

Sources proposed as contributers to the annihilation emission have included cosmic-ray 
interactions in the interstellar medium (Lingenfelter & Ramaty 1982), pulsars (Sturrock 1971) and 
(3 + decay products from radioactive nuclei (e.g., 56 Co, 44 Sc, and 26 Al) produced by supernovae, 
novae, or Wolf-Rayet stars (Clayton 1973; Ramaty & Lingenfelter 1979; Signore & Vedrenne 1988; 
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Woosley & Pinto 1988; Lingenfelter & Ramaty 1989). However, efforts to model the distribution 
of the 511 keV emission have been only moderately successful (see, e.g., Purcell et al. 1993). 
These models, consisting of components generally representing Galactic emission observed at 
other wavelengths (e.g., CO maps) or geometric configurations assuming that the emission of 
annihilation radiation follows a Galactic distribution of known sources, were not able to describe 
the OSSE maps adequately. It is evident that simple Galactic distributions are not sufficient 
to describe the most recent 511 keV data and that additional components are required. For 
example, the expected 511 keV line flux resulting from the radioactive decay of 26 Al represents 
only ~ 3 - 9% of the total Galactic 511 keV line flux (Purcell et al. 1997). 

In this paper, we explore the possibility that Sgr A East, a very prominent elongated radio 
structure surrounding (although slightly off-centered from) the Galactic nucleus, is a significant 
source of the annihilation radiation associated with the bulge feature. The motivation for doing 
so stems from the association of this radio object and the EGRET 7-ray source 2EG J1746-2852 
suggested by Melia et al. (1998). These authors show that the high-energy component of the 
EGRET spectrum can be accounted for by the decay of neutral pions (ir° — > 77) produced in p — p 
scatterings between ambient protons and a relativistic population accelerated by shocks within the 
Sgr A East shell. Since charged pions which decay leptonically (ir^ — > /r* 1 ^, with //^ — > e^feU^) 
are also produced, the rest of Sgr A East's spectrum, extending down to GHz energies, is filled 
uniquely by the bremsstrahlung, Compton and synchrotron interactions of the associated electrons 
and positrons. The VLA data, for example, are consistent with the synchrotron emissivity of the 
decay leptons within an equipartition magnetic field of « 10~ 5 G. 

If this scenario is correct, Sgr A East must be a site of copious positron production. Here, we 
investigate whether this process is sufficient to power the measured flux of annihilation radiation 
from the Galactic Center, and if not, whether theses observations suggest an evolutionary history 
of Sgr A East. 

The rest of the paper is organized as follows. Section 2 presents what is known about Sgr A 
East from observations as well as what we infer based on the results of Melia et al. (1998). This 
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section also includes a calculation of the ionization state of the gas in Sgr A East. The fate of 
positrons produced in Sgr A East is then considered in Section 3. Two different evolutionary cases 
for this environment are presented in Section 4. The characteristics of the annihilation radiation 
spectra for these cases are determined in Section 5 and compared to observations. Our concluding 
remarks are also given in Section 5. 

2. The Environment of Sgr A East 

2.1. The Basic Properties 

Sgr A East is an elliptical structure with nonthermal radio emission peaked at its periphery. 
Geometrically, the structure is elongated along the Galactic plane with a major axis of length 
10.5 pc and a center displaced from the apparent dynamical nucleus, Sgr A West, by 2.5 pc in 
projection toward negative Galactic latitudes. The actual distance between Sgr A West and the 
geometric center of Sgr A East has been estimated to be ~ 7 pc (Yusef-Zadeh & Morris 1987; 
Pedlar et al. 1989; Yusef-Zadeh et al. 1999). Although morphologically similar to supernova 
remnants (SNRs), the energetics (~ 4 x 10 52 ergs) of Sgr A East, based on the power required to 
carve out the radio synchrotron remnant within the surrounding dense molecular cloud, appear 
to be extreme compared to the total energy (~ 10 51 ergs) released in a typical supernova (SN) 
explosion. It has been suggested that Sgr A East may be the remnant of a tidally disrupted 
star (Khokhlov & Melia 1996). Regardless of what actually caused the initial explosion, recent 
observations of this region at 1720 MHz (the transition frequency of OH maser emission) have 
revealed the presence of several maser spots at the SE boundary of Sgr A East with a velocity 
of ss 50 km s -1 , and one near the Northern arm of Sgr A West at a velocity of 134 km s _1 
(Yusef-Zadeh et al. 1996). These observations are consistent with the presence of shocks produced 
at the interface between the expanding Sgr A East shell and the surrounding environment. The 
implied age of Sgr A East (r ~ 5pc/100 km s _1 ~5x 10 4 yr) is then consistent with those of 
typical SNRs. 
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Given its proximity to the Galactic Center, Sgr A East is bathed by intense IR and UV 
radiation fields associated with the central 1^2 parsecs of the galaxy (Telesco et al. 1998; 
Davidson et al. 1992; Becklin, Gatley & Werner 1982; Melia, Yusef-Zadeh & Fatuzzo 1998). 
Although there are no direct measurements of the properties of the ambient medium enveloped 
by Sgr A East, X-ray observations suggest an average electron ambient density of ~ 3 — 6 cm~ 3 
(Koyama et al. 1996; Sidoli & Mereghetti 1999). These results are consistent with an upper 
limit of ~ 30 cm -3 inferred by Melia et al. (1998). The temperature of the ambient medium 
is considerably less constrained. Observationally, the ISM within the inner 150 pc is comprised 
of several gas components with a range of temperatures. Inside the minicavity, the electron 
temperature lies in the range 4000 — 7000 K, and is comparable to the average value of 7000 K 
found for Sgr A West (Roberts, Yusef-Zadeh & Goss 1996). However, on a larger scale (i.e., tens 
of pc), the X-ray spectra observed with ASCA exhibit several emission lines from highly ionized 
elements that are characteristic of a ~ 10 keV thermal plasma. In particular, this region stands 
out in the intense 6.7 keV Ka-transition of He-like Fe (Yamauchi, et al. 1990; Koyama, et al. 
1996). These spectra indicate a common origin for the plasma in the entire region, rather than 
a superposition of individual sources. Assuming an expansion velocity comparable to the sound 
speed for kT ~ 10 keV, Koyama et al. (1996) infer an expansion age of ~ 10 5 yr, which is similar to 
that of Sgr A East. Thus, both the X-ray plasma and the relativistic synchrotron source may have 
been produced by the same explosive event ~ 0.5 — 1 x 10 5 yr ago. The BeppoSAX observations of 
the Sgr A Complex suggest somewhat different characteristics, in which the high-energy data are 
well fit by the sum of two thermal models with kT ~ 0.6 and 8 keV, with the lower temperature 
plasma well correlated to Sgr A East (Sidoli & Mereghetti 1999). 

2.2. Other Inferred Characteristics 

Adopting the scenario of Melia et al. (1998) invoked to account for the Galactic Center 
EGRET observations, we assume that pion-production occurs within the Sgr A East shell as a 
result of collisions between shock-accelerated and ambient protons (see also Markoff, Melia & 
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Sarcevic 1997). A crucial aspect of this mechanism is the decay of charged pions to muons, and 
subsequently, to relativistic electrons and positrons. The rate at which these leptons are produced 
is linked directly to the rate at which neutral pions are produced (since the charged pions and 
neutral pions are both byproducts of the same scattering events), and is therefore well constrained 
by the EGRET observations (Markoff, Melia & Sarcevic 1999). In addition, since the relativistic 
leptons radiate synchrotron emission, their energy distribution is well constrained by the radio 
properties of Sgr A East. These observations also provide an indirect measure of the magnetic 
field strength in Sgr A East (see below). 

In order to keep the analysis as concise as possible while at the same time capturing the 
important characteristics of the Sgr A East environment, we assume a uniform and homogeneous 
shell geometry, with an inner radius of 4 pc and an outer radius of 5pc, and a centroid about R = 
7 pc from the Galactic Center. The IR and UV radiation fields bathing this region are assumed 
to have blackbody spectra with temperatures Tj r = 100 K and T uv = 30, 000 K, respectively, and 
energy densities 

- - - 44 x 10 "'° ^ ^ {ttwrd (£) " 2 ■ (1) 

where L is either Lj r = 10 7 L Q , or L uv = 2 x 10 7 L Q . Furthermore, we will assume that Sgr A 
East has an age of 75, 000 years and is comprised mostly of hydrogen with an average density 
riH = 5 cm -3 , and helium with an average density riHe = 0.5cm -3 . The temperature profile T 
within the Sgr A Complex must reflect the superposition of the various components discussed 
above. For our purposes here, however, we assume that the ambient plasma where the positrons 
are produced may be ascribed a single characteristic value of T at any given time. The positrons 
are assumed to be produced within the Sgr A East shell in accordance with Case 2 of Melia et al. 
(1998), the one that produced the best fit to the gamma-ray data. The injection function I + for 
these particles is shown in Figure 1 as a function of the positron energy E + . The total production 
rate, found by integrating this function over all values of E + , is equal to 2.7 x 10~ 18 cm -3 s _1 . 
The corresponding magnetic field strength, inferred from the VLA spectrum, is B = 1.1 x 10 -5 G, 
which is within a factor of 2 of its equipartition value. 
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One of the processes by which annihilation occurs is charge exchange between positrons and 
atoms. Our analysis of the electron-positron annihilation considered below will therefore require 
an assessment of the ionization fractions in the Sgr A East shell, which are set by the balance 
of total recombination and ionization rates for hydrogen and helium. For simplicity, we assume 
that all atoms are ionized from their ground state. We will confirm a posteriori that the ambient 
plasma in Sgr A East is highly ionized and that therefore the approximations we make here do 
not greatly affect the results. We can therefore ignore the trace presence of neutral helium as the 
charge exchange process is dominated by the more prevelant singly ionized species of this atom. 
In addition, the number density of ambient electrons is well-approximated by the constant value 
n e = 6 cm -3 . The recombination rates are therefore independent of the ionization fractions of 
the species comprising the plasma, and hence the ionization equations for hydrogen and helium 
decouple and reduce to the following: 

n 



n HII R HI + C HI > 
n HeII ^Rff 111 



n HeIII R HeII + fjHell > 



(2) 



(3) 



where Rfi is the radiative recombination rate to state i and R\f and C\f are the radiative and 
collisional ionization rates from the ground state. Collisional recombination can be ignored for the 
conditions found in Sgr A East. 

Since the electron self-collision rate 



(Spitzer 1956) is much larger than the radiative recombination rates in the Sgr A East environment 
(see Figure 2), the free electrons thermalize prior to recombining. We can therefore use the 
collisional ionization rates for HI and Hell in terms of the ambient electron temperature T and 
the corresponding ground-state ionization energy x from Bell et al. (1983). These represent the 
best fits to the observation and theory, and take the form 

5 

C lf = n e exp(- X /kT) (kT/ X ) 1/2 £ a n [log(kT/ X )] n , (5) 

n=0 



when x/10 < kT < 10 X , and 

C lf = n e (kT/ X )- 1/2 



7 HkT/ X ) + J2Pn(x/kTy 

n=0 



(6) 



when kT > 10%. The values of the coefficients are given in Tables 1 and 2. 

Table 1 
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HI 


2.3742E-08 


-3.6866E-09 


-1.0366E-08 


-3.8010E-09 


3.4159E-09 


1.6834E-09 


Hell 


3.4356E-09 


-1.6865E-09 


-6.9236E-10 


9.7863E-11 


1.5591E-10 


6.2236E-11 
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HI 


2.4617E-08 


9.5986E-08 


-9.2463E-07 


3.9973E-06 








Hell 


3.0772E-09 


1.1902E-08 


-1.1514E-07 


5.0489E-07 





The radiative ionization rate from the ground state is given by the expression 

roo 

Rif = / cr 1 f(e)cn t de , 
where n t is the differential number density of ionizing UV photons, normalized so that 



(7) 



(8) 



The ionization cross-section for hydrogenic atoms (e.g., HI and Hell) in their ground state is given 
by 



<^i/(e) 



8h 3 



X 

S 3 / 2 TT 2 e 2 m^cZ 2 \ e 



'' ' : 1 _ e -2n(a/f3) 



(9) 



where Z is the atomic number, a is the fine structure constant, (3 = v/c is the ratio of the free 
electron speed to that of light, and the term in square brackets is an analytical expression for the 
gaunt factor (Spitzer 1978). 



The total radiative recombination rate for hydrogenic atoms is 



(10) 
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where v and E are, respectively, the speed and kinetic energy of the electrons prior to capture, 

2 4 he 2 x 2 9 

is the cross-section for radiative capture to the i-th. state, and the brackets <> denote an average 
over the thermal electron distribution function. The sum over i may be truncated at i = 20 
without introducing any appreciable error. For simplicity, the gaunt factors g are set equal to 
unity, therefore introducing a possible error that ranges from ~ 3% at T = 10 4 K to ~ 20% at 
temperatures in excess of 10 6 K. We note, however, that our final results are not subject to the 
large errors incurred at temperatures > 10 6 K. 

Values of the rate coefficients as a function of the ambient temperature T are shown in Figure 
2. Note that the radiative recombination rates are lower at all temperatures than the electron 
self-collision rate, as discussed above. For these conditions then, the ionization equations are easily 
solved to determine the atomic number densities. The results of these calculations are shown in 
Figure 3, and these help to justify the assumption stated earlier that the plasma is highly ionized 
for the environment in the Sgr A East shell. 



3. The Positron Population 

In this section, we calculate the positron distribution, and determine whether these particles 
thermalize before annihilating. We then consider how to calculate the present positron number 
density for a given evolution of the Sgr A East environment. 

3.1. The Positron Annihilation Rates 

Positron annihilation in a hydrogen and helium plasma can occur via three dominant 
mechanisms: direct (in-flight) annihilation; radiative capture; and charge exchange. The first of 
these processes, in which a free positron annihilates with an ambient electron without first forming 
positronium, results in the production of two photons. The direct annihilation rate for a positron 
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of kinetic energy E + interacting with an isotropic Maxwell-Boltzman distribution n e f m b(E_) of 
electrons (where E_ is the electron kinetic energy) is given by the expression 

- 1 dfi 



R da {E+) = cn e J f mb (E-)dE- J Jlf kin (E + , E-) a da (E+,E-) , 



where 



1/2 



(12) 



(13) 



fkin = [pi ai0i (i - /i 2 ) - w-P+v 

accounts for the kinematics of the two-particle motion, and 

[ 2E ± m e c 2 + El 1 1/2 
;± [(mc 2 ) 2 + 2E ± m e c 2 + El\ 1 ' 

are the positron/electron speeds in units of c. The cross-section as a function of the relative 
Lorentz factor 7 = (1 — /3 2 ) -1 can be determined using the plane-wave approximation and is found 
to be 



&da = 

where 



2tt 



7 2 + 4 7 + 1 



r 



7 + 1 V 

(l-^_M) 2 -(l-^)(l-/9l) 



1/2 



(15) 



(16) 



(l-/3 + /3_/x) 2 

is the relative speed (in units of c) between the electrons and positrons (Coppi & Blandford 1990), 
ro is the classical electron radius, and the term in the first set of square brackets takes into account 
the Coulomb attraction. 

Radiative capture produces positronium in either a singlet spin-0 state (0.25 probability) or 
a triplet spin-1 state (0.75 probability). Annihilation from the singlet state produces two line 
photons whereas annihilation from the triplet state produces three continuum photons. The 
overall capture rate is given by summing over the radiative capture rates into each state i: 

n dn 



R rc (E + )=Y J cn e J f mb (E_)dE_ J^f km (E + ,E_) a rc (E + ,E_) , 



(17) 



where the cross-section is easily generalized from Equation (11) by replacing m e with the reduced 
mass m e /2 and noting that the ionization energy from the ground state of positronium is % = 6.8 
eV. 
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The presence of Hell implies that charge exhange (e + + Hell — > Ps + Helll) is an important 
mechanism by which positronium is formed. The rate at which this process occurs is given by 

R ce {E + ) = n H eii (3+ca ce (E + ) . (18) 

To our knowledge, no generally accepted theoretical or experimental value exists for the 
cross-section of this process. However, based on measurements of the cross-sections for the e + + 
HI ^ Ps + HII (Sperber et al. 1992) and the e+ + Hel -> Ps + Hell reactions (Overton et al. 
1993), it is clear that the desired cross-section is peaked at 54.4 eV (the ground-state ionization 
energy for Hell) with a value there of ~ ttoq. Assuming an energy profile similar to that of the e + 
+ Hel — > Ps + Hell cross-section, we may therefore adopt an empirically based cross-section of 
the form 

a ce (E + ) = 2.4 x 1(T 16 cm 2 (|± - l) e^~^) , (19) 
for positron kinetic energies E + > E\ = 28 eV. 

The annihilation rates are shown as functions of the positron kinetic energy in Figures 4-6 
for three different values of the ambient temperature. The dotted lines, which represent the 
Maxwell-Boltzman functions / m j,(i?_) at each specified temperature, will be introduced in § 4. 



3.2. The Positron Cooling Rates 



We present here the cooling rates for a positron moving through Sgr A East with a Lorentz 
factor 7+ = [1 — /3 2 .] -1 . The dominant energy loss mechanisms are: (1) Synchrotron losses at a 
rate 



dE s 
dt 

(2) Compton scattering losses at a rate 



4cj tc 2 2 R 2 . 



dE n 
~dt 



m e c 2 + m e c 2 + j + kTi 



where x uv = 2.7^ + kT uv / m e c 2 , Xi r = 2. 77+ /cTj r /m e c 2 , and 



a c (x) 



-<jt 



^ ( 2X{1+X) - HI + 2,)) + -1 ln(l + 2x) - , 1 + 3X 
x 3 V 1 + 2x V 7 2x V ' (1 + 2xf 



(20) 



(21) 



(22) 
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(3) Bremmstrahlung emission losses at a rate 



(JL L! -4ar 2 m e c 3 (2n H + 6n He ) ln(2 7+ - l/3) 7+ ; (23) 



dt 

and (4) Coulomb losses at a rate of 

dE c 3 f £ /7T1 _^ J771 _ f 1 dfi 

dt 

where 

3ctt(7- - 7+) 



m e c 3 n e J f mb (E )dE j f kin < aAj > , (24) 



< 0-A7 >= 



64e 2 



(4 _ A + 4 ,n 2) (t*£ V- 2 (?^LL) + - S 4 £ + £ + 



p2 1 l p2 £ 2 I £ 4 3 e 4 e 4 

(25) 

is given in terms of the parameters e 2 = (1 +j)/2 and p 2 = e 2 — 1, where 7 is the relative Lorentz 
factor (Coppi and Blandford 1990). The corresponding cooling rates, defined by the general 
expression R = E^ 1 (dE/dt), are plotted in Figure 7 for an ambient temperature of 10 5 K. We note 
that only the Coulomb cooling rate depends on the temperature, with the low-energy turnover 
occur ing at an energy of ~ kT. 



3.3. The Fate of Positrons in Sgr A East 

It is clear from the results shown in Figures 4-7 that the total positron cooling rate is much 
greater than the total positron annihilation rate over the entire energy range. In addition, for 
energies below ~ 10 8 eV, the total cooling rate increases with decreasing particle energy. Taken 
together, these results imply that particles produced in Sgr A East with energies < 10 8 eV 
thermalize to the same temperature as the ambient electrons on a time scale equal to the inverse 
of their initial cooling rate, and do so before they annihilate. These results are summarized in 
Figure 8, which shows the total annihilation (short-dash) and cooling (long-dash) rates assuming 
an ambient temperature of T = 10 5 K along with the injection function (solid line, right-hand 
scale) for the Sgr A East positrons. The two horizontal dotted lines mark the inverse of 10 4 
and 10 5 years, and the energy at which they first intersect the injection function then yields the 
threshold energy below which particles (i.e., those to the left of the corresponding vertical lines) 
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have had sufficient time to thermalize within that time scale. A plot of the threshold energy E t 
as a function of time ranging between 10 4 and 10 5 years is shown in Figure 9. Figure 10 then 
illustrates the rate per unit volume with which positrons thermalize in Sgr A East (found by 
integrating the injection function up to the threshold energy) as a function of its age. This value 
falls well short of the total rate per unit volume with which positrons are produced in Sgr A East 
(2.7 x 1CT 18 cm~ 3 s _1 ) since E t is smaller than the energy at the peak of the injection function. 

In other words, Sgr A East is apparently producing positrons at a rate that is currently much 
higher than that with which they cool and thermalize. The positrons in this remnant are being 
"stockpiled." 

4. Evolutionary Models of Sgr A East 

The number density of thermal positrons populating Sgr A East as a function of its age t* is 
found by integrating the evolution equation 

- [n+ (< R da > + < R ce > + < Rrc >)] , (26) 

where I + is the injection function shown in Figure 1 and the terms in the <> brackets represent 
the annihilation rates (da - direct annihilation; ce - charge exchange; rc - radiative capture) 
averaged over the thermal positron distribution function. Note that the term in the first square 
brackets represents the rate at which positrons are thermalized in Sgr A East. 

We first consider a static evolutionary model for Sgr A East where both the temperature 
and thermalization rate are constant in time. For this case, the number density of positrons 
increases with age until the annihilation rate matches the thermalization rate. Assuming that 
all thermal positrons form positronium before annihilating, each positron produces an average 
of 0.25 x 2 + 0.75 x 3 = 2.75 annihilation photons. A reasonably good estimate for the ensuing 
flux at Earth, found by multiplying the product of the injection rate (which is clearly larger than 
the thermalization/annihilation rate) and the Sgr A East volume by 2.75 to find the total rate 
of annihilation radiation photons produced in Sgr A East, and then dividing by 4irD 2 (where 



dn 



+ 



dt 



E t {t*) 



I+(E)dE 
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D = 8.5 kpc is the distance to the Galactic Center), yields a value of 6.5 x 10 6 cm 2 s 1 . This 
number is clearly much smaller than the measured flux arriving at earth from the bulge. 

Is it possible then that a time-dependent positron production rate in Sgr A East could 
account for the current high rate of annihilation? Let us consider the possibility that positrons 
were stockpiled in Sgr A East during its infancy, and are now being depleted at a rate exceeding 
the current thermalization rate. This feat can be accomplished in two ways. First, the positron 
injection rate in Sgr A East may have been much higher earlier in its life. Indeed, the energy 
content of the relativistic proton population for Case 1 of Melia et al. (1989) is 7.7 x 10 51 ergs, or 
roughly 10% of the energy associated with Sgr A East. (We note, however, that this value is really 
an upper limit, as it was derived assuming steady state conditions.) In that case, the number 
density of positrons populating Sgr A East would now be ~ 10~ 2 cm~ 3 . In contrast, for the static 
evolutionary model discussed above, the number density of positrons injected into Sgr A East 
during the assumed lifetime of 75,000 years is ~ 6 x 10~ 6 cm -3 . This difference suggests that the 
initial injection rate may have been as high as 10 4 times greater than what is observed today. 

We therefore consider two cases for the evolution of the injection rate, both of the form 

/+(*) = A /+(**) exp(-t/r), (27) 

where A = 10 2 for Case A and A = 10 4 for Case B. The value for r is then well-defined for the 
assumed age t* of 75,000 years. 

Another reason why positrons may have been stockpiled in Sgr A East is that its temperature 
decreases with age. This can be seen by considering the annihilation rates averaged over 
the thermal positron distribution as a function of temperature. As shown in Figure 11, the 
averaged charge-exchange rate is strongly peaked at ~ 10 5 K, where the overlap between the 
Maxwell-Boltzman distribution function and the charge-exchange rate produces a clear maximum 
(see Figures 4-6). It is intriguing that the lower limit of 0.94 for the positronium fraction 
observed from the Galactic Center indicates that the temperature where the annihilation radiation 
is produced must lie in the range between 3.2 x 10 4 K and 3.2 x 10 5 K. This result is clearly seen 
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from the plot of the positronium fraction 

f Ps = < ^ce > + < Rrc > (2g) 

< R ce > + < Rrc > + < Rda > 

versus temperature, as shown in Figure 12. These results are somewhat consistent with the 
observed line-width of 2.5 keV, since the expected relation 

T * 3 - 2 * 10 " {wkT • < 29 > 

(Wallyn et al. 1996) then implies a temperature of ~ 2.8 x 10 4 K. We caution, however, that the 
above line- width relation is somewhat suspect (as noted by Wallyn et al.). We therefore assume 
that the initial temperature of Sgr A East may have been 10 7 K and that it decays exponentially 
so that the temperature will be 10 5 K at an age of 10 5 years. 



5. Results and Discusson 

We have integrated Equation (26) for the two evolutionary cases (A and B) discussed above. 
The number densites as a function of age for these situations are shown in Figure 13 as solid lines. 
The dotted lines represent the number densities stockpiled in the event that the annihilation rates 
are zero. To gain insight into these results, we plot the rate per unit volume at which positrons 
are thermalized (dotted line) and annihilated (solid and dashed lines) as a function of age for case 
B in Figure 14. It is clear that the positrons are indeed stockpiled when the thermalization rates 
are high, and this stockpile is depleted when the charge-exchange annihilation rate peaks. 

The ensuing annihation spectrum is easily determined by noting that the observed flux of line 
photons is given by 

Fi = (2 < Rda > +0.5 < R rc > +0.5 < R ce >) , (30) 

and the observed flux of continuum photons is given by 

F c = (2.25 < R rc > +2.25 < R ce >) , (31) 
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where N + is the total number of thermal positrons in Sgr A East. The total flux observed as a 
function of age is presented in Figure 15 for the two cases discussed above. As expected, the flux 
peaks at around 75, 000 years, where the positrons are just beginning to be depleted. 

It is clear by comparing the values in Figure 15 to the observed flux from the bulge that 
even ideal evolutionary conditions cannot make up for the fact that only a small fraction of the 
injected positrons will have thermalized and are annihilating. This point is emphasized by the 
requirement of a very efficient cooling mechanism for the conditions in Sgr A East to account for 
the observed annihilation radiation flux, as may be seen by making the ad hoc assumption that all 
of the injected positrons thermalize. In that (extreme) case, the resulting flux curve (as shown in 
Figure 16 for the case with A = 10 2 ) would agree well with the observations. 

Our first conclusion is therefore that Sgr A East is not (by itself) the source of annihilation 
radiation from the Galactic Center, unless some efficient cooling mechanism (not considered here) 
is thermalizing the positrons on a time scale shorter than the age of this remnant. 

However, since Sgr A East's age is much shorter than the age of the Galaxy, it is reasonable 
to ask whether the process that produced this relatively unique remnant was a singular event. 
If not, then the accumulated "stockpiles" of positrons may now be at the required level to 
account for the observed annihilation rate. Unfortunately, there is very little to go by in terms of 
phenomenologically pinning down the rate at which these explosive events have occurred. But we 
can at least use the properties of the bulge component of the e + e~ annihilation radiation to infer 
what the minimum explosion rate must have been. 

Using the required positron density, and assuming that no special cooling rate, beyond what 
we have considered here, is acting on the positrons, we estimate that ~ 30 explosions (like the one 
that produced Sgr A East) must have occurred over the past one million years if A = 10 2 . The 
number of events drops as the value of A increases, reaching about one event per million years 
when A = 10 4 . This is interesting in view of the fact that Sgr A East's age is estimated to be 
about 75,000 years, which means that indeed we may be witnessing the most recent of a series of 
catastrophic events that occur on average about once every 30,000 to 1 million years. This lends 



-18- 



indirect support to the idea (Khokhlov and Melia 1996) that Sgr A East may have been produced 
by the tidal disruption of a ~ 2 M Q star that ventured within 10 Schwarzschild radii or so of the 
central black hole. Given the stellar density in the vicinity of Sgr A*, such events are expected to 
occur roughly once every 10,000 to 100,000 years. The bulge component of the e + e~ annihilation 
radiation may therefore be the best evidence we have right now for the existence of accumulated 
relics of this ongoing stellar disruption process. 
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Fig. 1. — The injection function I + (E + ) for positrons produced via charged pion decay in Sgr A 
East, assuming conditions adopted in Case 2 of Melia et al. (1998). The total rate per unit volume 
at which positrons are produced, determined by integrating this function over the positron kinetic 
energy, is equal to 2.7 x 10~ 18 cm" 3 s -1 . 
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Fig. 2. — Recombination and ionization rates for hydrogen and helium in Sgr A East as a function 
of the ambient electron temperature. 
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Fig. 3. — The number density of neutral and ionized hydrogen and of singly and doubly ionized 
helium in Sgr A East as a function of the ambient electron temperature. 
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Fig. 4. — Annihilation rates as a function of positron kinetic energy for an ambient temperature of 
10 4 K. 
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Fig. 5. — Annihilation rates as a function of positron kinetic energy for an ambient temperature of 
10 5 K. 
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Fig. 6. — Annihilation rates as a function of positron kinetic energy for an ambient temperature of 
10 6 K. 



- 28 - 



n i i r 



-2 



-4 



-6 



K 

S 

o 
o 



o2 -10 

o 



10 5 K 



-12 



-14 



-16 



Synchrotron 
Compton 
Bremmstrahlung 
Coulomb 




2 4 6 8 

log 10 Positron Kinetic Energy (eV) 



Fig. 7. — Cooling rates as a function of positron kinetic energy for an ambient temperature of 10 5 
K. 
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Fig. 8. — A direct comparison between the total annihilation and cooling rates, and the injection 
function in terms of the positron kinetic energy for an ambient temperature of 10 5 K. 
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Fig. 9. — The threshold energy below which positrons thermalize on a time scale shorter than the 
age of Sgr A East as a function of the remnant's age. 
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10. 



The thermalization rate as a function of age for Sgr A East. 
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Fig. 11. — The positron annihilation rates averaged over the positron distribution as a function of 
temperature. 
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Fig. 12. — The positronium fraction as a function of temperaure. 
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Fig. 13.— 



The number density of thermal positrons for the two cases considered. 
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Fig. 14. — A comparison between the rates per unit volume of thermalization and annihilation as 
a function of age for A = 10 4 . 
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Fig. 15. — The flux (as observed at Earth) of the total number of annihilation photons for both 
A = 10 2 and A = 10 4 as a function of the age of Sgr A East. The flux peaks shortly after the 
annihilation rate due to charge exchange exceeds the thermalization rate (see Figure 14). 
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Fig. 16. — The flux as observed on Earth of the total number of annihilation photons for A = 10 2 
under the ad hoc assumption that all injected positrons thermalize on a time scale that is much 
shorter than the age of Sgr A East. 



